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A mapping theorem leading to exact many-body dynam- 
ics of impenetrable bosons in one dimension reveals dark and 
gray soliton-like structures in a toroidal trap which is phase- 
imprinted. On long time scales revivals appear that are be- 
yond the usual mean-field theory. 
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Dark and gray solitons are a generic feature of the 
nonlinear Schrodinger equation with repulsive interac- 
tions, and several calculations of their dynamics based on 
the Gross-Pitaevskii (GP) equation have appeared 
as well as experiments demonstrating their existence in 
atomic BECs Since the GP equation is a nonlin- 

ear approximation to the more exact linear many-body 
Schrodinger equation, this raises the question of how ob- 
served solitonic behavior arises in a theory which is linear 
at a fundamental level. Here this issue will be examined 
with the aid of exact many-body solutions. It has been 
shown by Olshanii |^ that at sufficiently low tempera- 
tures, densities, and large positive scattering length, a 
BEC in a thin atom waveguide has dynamics which ap- 
proach those of a one-dimensional (ID) gas of impene- 
trable point bosons. This is a model for which the exact 
many-body energy eigensolutions were found in 1960 us- 
ing an exact mapping from the Hilbert space of energy 
eigenstates of an ideal gas of fictitious spinless fermions 
to that of many-body eigenstates of impenetrable, and 
therefore strongly interacting, bosons pO| , pd[ |. The term 
"Bose-Einstein condensation" is used here in a general- 
ized sense; it was shown by Lenard and by Yang and 
Yang that for the many-boson ground state of this 
system, the occupation of the lowest single-particle state 
is of order y/N where N is the total number of atoms, in 
contrast to N for usual BEC. Nevertheless, since ^ 1 
and the momentum distribution has a sharp peak in the 
neighborhood of zero momentum this system shows 
strong coherence effects typical of BEC. The response of 
a BEC of this type to application of a delta-pulsed optical 
lattice was recently calculated by Rojo et al. jl^, using 
the Fermi-Bose mapping theorem as an exactly 

calculable model of dynamical optical lattice behavior. 
They found spatial focussing and periodic self-imaging 
(Talbot effect), which decay as a result of interactions. 
This decay is absent in the GP approximation and there- 
fore serves as a signature of many-body interaction effects 
omitted in GP. 



In this Letter we examine the appearence of dark 
soliton-like structures using the model of a ID hard-core 
Bose gas in a toroidal trap, or ring, with cross section so 
small that motion is essentially circumferencial [p"5[-p^. 
The Fermi-Bose mapping is employed to generate exact 
solutions for this problem. We identify stationary solu- 
tions which reflect some properties of dark solitons from 
the GP theory when the ring is pierced at a point by an 
intense blue-detuned laser. We also present dynamical 
solutions when half of an initially homogeneous ring BEC 
is phase-imprinted via the light-shift potential of an ap- 
plied laser, leading to gray soliton-like structures whose 
velocity depends on the imposed phase-shift [HJ|]. Such 
structures are apparent for times less than the echo time 
Te = L/c, with L the ring circumference and c the speed 
of sound in the BEC. On longer time scales the dynamics 
becomes very complex showing Talbot recurrences which 
are beyond the GP theory. 

Time- dependent Fermi-Bose mapping theorem: The 
original proof was restricted to energy eigenstates, 

but the generalization to the time-dependent case is al- 
most trivial. The Schrodinger Hamiltonian is assumed to 
have the structure 
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where Xj is the one-dimensional position of the jth par- 
ticle and V is symmetric (invariant) under permutations 
of the particles. The two-particle interaction potential is 
assumed to contain a hard core of ID diameter a. This 
is conveniently treated as a constraint on allowed wave 
functions 'ijj{xi, ■ ■ ■ , ; t) : 



i/- = if |a 



Xk\ < a 



I < j < k < N, (2) 



rather than as an infinite contribution to V, which 
then consists of all other (finite) interactions and ex- 
ternal potentials. The time-dependent version starts 
from fermionic solutions ippixi, ■ ■ ■ , XN;t) of the time- 
dependent many-body Schrodinger equation (TDMBSE) 
Hip = ihdip/dt which are antisymmetric under all parti- 
cle pair exchanges Xj ^ Xk, hence all permutations. As 
in the original theorem [p^ , define a "unit antisymmetric 
function" 

^(a;i, • • • ,a;Ar) = Y]_ sgn{xk-Xj), (3) 

i<j<k<N 
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where sgn(a;) is the algebraic sign of the coordinate dif- 
ference it is +1(-1) if a; > 0{x < 0). For 
given antisymmetric ipp, define a bosonic wave function 
i/jB by 

iPb{xi, - ■ ■ ,XN;t) ^ A{xi, - ■ ■ ,XN)lpFixi, - ■ ■ ,XN;t) (4) 

which defines the Fermi-Bose mapping. ipB satisfies the 
hard core constraint (2) if ipp does, is totally symmetric 
(bosonic) under permutations, obeys the same bound- 
ary conditions as ipp, e.g. periodic boundary condi- 
tions on a ring, and HipB — ihdipB/dt follows from 

Exact solutions for impenetrable point bosons: The 
mapping theorem leads to explicit expressions for all 
many-body energy eigenstates and eigenvalues of a ID 
scalar condensate (bosons all of the same spin) under 
the assumption that the only two-particle interaction 
is a zero-range hard core repulsion, represented by the 
a limit of the hard-core constraint. Such so- 
lutions were obtained in Sec. 3 of the original work 
for periodic boundary conditions and no exter- 
nal potential. The exact many body ground state 
was found to be a pair product of Bijl-Jastrow form: 
■^o = const, fl^^^ I sin[7rL~^(a;i — Xj)]\. In spite of the 
very long range of the individual pair correlation fac- 
tors I sin[7rL~^(a;i — Xj)]\, the pair distribution function 
D{xij), the integral of |'0oP over all but two coordinates, 
was found to be of short range: D{xij) = 1 — j'^{-Kpxij), 
with jo(C) = sin^/^, the spherical Bessel function of or- 
der zero. The system was found to support propagation 
of sound with speed c = nfip/m where p ~ N/L, the ID 
atom number density. 

To generalize to the time-dependent case, assuming 
that the many-body potential of Eq. (1) is a sum of 
one-body external potentials V{xj,t), one generalizes 
the time- independent determinantal many-fermion wave- 
function |10i to a determinant 



i>F{xi, 



N 



■ ,XN',t) — C det (j)i{xj,t), 



(5) 



of solutions <j)i{x,t) of the one-body TDSE in the exter- 
nal potential V{x, t). It then follows that t/jj? satisfies the 
TDMBSE, and it satisfies the impenetrability constraint 
(vanishing when any Xj = xi) trivially due to antisymme- 
try. Then by the mapping theorem -05 of Eq.(4) satisfies 
the same TDMBSE. 

Dark solitons on a ring: Consider N bosons in a tight 
toroidal trap, and denote their ID positions measured 
around the circumference by Xj . This is equivalent to the 
exactly-solved model of N impenetrable point bosons 
in ID with wave functions satisfying periodic boundary 
conditions with period L equal to the torus circumfer- 
ence, and the fundamental periodicity cell may be cho- 
sen as —L/2 < Xj < L/2. However, the rotationally 
invariant quantum states of this problem do not reveal 



any dark soliton-like structures. To proceed we therefore 
consider the case that a blue-detuned laser field pierces 
the ring at a; = by virtue of the associated repulsive 
dipole force: The light sheet then provides a reference 
position for the null of the dark soliton. Assume that the 
light sheet is so intense and narrow that it may be re- 
placed by a constraint that the many-body wave function 
(hence the orbitals (pi) must vanish whenever any Xj — 0. 
Then the appropriate orbitals (j)i{x) are free-particle en- 
ergy eigenstates vanishing at a; = and periodic with 
period L. The complete orthonormal set of even-parity 
eigenstates (f)^^^ and odd-parity eigenstates (piT^ are 

(pl^\x) = v/27Isin[(2n- l)7r|a;|/L], 

(t)i^\x) = y/2/Lsm{2m:x/L), (6) 

with n running from 1 to c». The odd eigenstates are the 
same as those of free particles with no x = constraint, 
since these already vanish at a; = 0. However, the even 
ones are strongly affected by the constraint, their cusp 
at a; = being a result of the impenetrable light sheet at 
that point. If one bends a ID box —L/2 < x < L/2 with 
impenetrable walls into a ring, identifying the walls at 
±L/2, then those particle-in-a-box eigenfunctions which 
are even about the box center become identical with the 
(f>l^\ and their cusp results from the nonzero slope of 
these functions at the walls. The A''-fermion ground state 
is obtained by inserting the lowest N orbitals (6) into the 
determinant (5) (filled Fermi sea). Assume that N is odd. 



Since 



is lower than 



this Fermi sea consists of 



the first {N + l)/2 of the </>l+^ and the first {N - I)/2 
of the (j}n . The iV-boson ground state is then given by 
(4). Since — 1, its one-particle density p{x) is the 
same as that of the A^-fermion ground state, the sum of 
partial densities contributed by all one-particle states in 
the Fermi sea. Thus it is the sum of 



p^+^x) = 



N + 1 sin[2(A^+ f)7ra;/L] 
2L 2Lsin(27ra;/L) 



(7) 



and 

p^-\x) 



N - 1 sin[(A - 1)ttx/L] cos[{N - 3)nx/L] 



2L 



Lsm{2Trx/ L) 



(8) 



In the thermodynamic limit N ^ oo, L ^ oo, N/L p 
for fixed x, p*-^^ each contribute half of the total density 
p{x): 



p{x) ^ p[l -jo{2TTpx)]. 



(9) 



Since jo(0) — 1, p{x) vanishes at a: = and approaches 
the mean density p over a healing length Lh — l/2p with 
damped spatial oscillations about its limiting value. This 
differs in detail from the density poo tanh^(a;/w) of a GP 
dark soliton [Ert, with poo the background density and 
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w the corresponding healing length, but has some quali- 
tative similarity. However, it is only the odd component 
p^-^\x) w p{x)/2 which has the feature of a dark soli- 
ton that the corresponding odd orbitals have a tt phase- 
jump at x = (and also at a; = ±L/2 to obey the peri- 
odic boundary conditions). But the odd and even com- 
ponents can never be separated physically, so the odd 
dark soliton-like component is always accompanied by 
the even non-soliton component. 

Next, suppose that the light-sheet is turned off at t = 
by removing the constraint that the wave function van- 
ish at x = 0. The solution of the TDMBSB for the 
many-boson system is then given by (4) where the Slater 
determinant (5) is built from the first {N + l)/2 of the 
(l)'h^\x,t) and the first (iV- l)/2 of the where 
these time-dependent orbitals are solutions of the single- 
free-particle TDSE which (a) reduce to the orbitals (6) at 
t = 0, and (b) satisfy periodic boundary conditions with 
periodicity cell —L/2 < x < L/2. The odd solutions are 
trivial: Since these never "see" the x = constraint even 
for t < 0, they differ from the odd orbitals (6) only by 



time-dependent phase shifts: 4'ri~\x,t) — (j)';^'{x)e 
with uj„ = hk^/2m and fc„ = 2mr/L. It follows that 
p^~\x,t) is time- independent, and given in the thermo- 
dynamic limit by 



p(-)(x,t)^(p/2)[l-jo(27rpx)]. 



(10) 



This further reinforces our view that the odd compo- 
nent of the density shares features of a dark soliton. The 
even-parity orbitals t/)!^'' {x, t) are complicated since the 
removal of the light sheet constitutes a large, sudden 
perturbation. Indeed, the periodic even-parity solutions 
of the free-particle Schrodinger equation are Xp^\^) ^ 
\J[2 — Spo)/L cos{2pTTx/L) with p = 0, 1, 2, • • •, and these 

are very different from the solutions <j)i^\x) with the 
X = constraint [Eq. (6)]. Nevertheless, since the 
Xp^\x) are complete for the subspace of even-parity, spa- 
tially periodic functions, one can expand the (f>i^^ {x, t) in 
terms of the Xp^'^ i^) i which evolve with time-dependent 
phases e~*'^p* with cUp — hkp/2m and kp = 2pTi/L. One 
finds 



2(2n - 1) PI ^{2- 5po) cos{kpx)e- 



L ^ 

p=0 



(2n- 1) 



(11) 



p(+^(a;,t) is the sum of absolute squares of the first 
(TV -I- l)/2 of the sums (11), generalizing (7). Adding 
the time-independent expression p^~^(x,t), given in the 
thermodynamic limit by (10) or exactly by (8), one finds 
the time-dependent total density p{x,t). There are two 
important time scales: One is the Poincare recurrence 
time Tr. Noting that Up in (11) is proportional to p^, one 



finds that all terms in the sum are time-periodic with pe- 
riod Tr = toL^/ttTi, which is therefore the recurrence time 
for the density and in fact all properties of our model . 
The other important time is the echo time t^, the time 
for sound to make one circuit around the torus. Recalling 
that the speed of sound in this system is c = nhp/m [ p^ , 
one finds Tg = Tr/N. For t << Tg after the constraint 
is removed, the initial density develops sound waves that 
propagate around the ring, and that we examine below 
in the context of phase-imprinting. For t > the evolu- 
tion is very complex, but complete recurrences occur for 
times t = riTr with fractional revivals in between. 

Gray soliton formation by phase-imprinting: Consider 
next a toroidal BEG in its ground state to which a phase- 
imprinting laser is applied over half the ring at i = 0. 
This is described by the single-particle Hamiltonian 



N 



-^^-nAosit)six,) 



(12) 



where S{x) — d{L/A — \x\), i.e., unity for — i/4 < a; < 
L/4 and zero elsewhere. This is the technique used in re- 
cent experiments [^,^ , here idealized to a delta- function 
in time and to sharp spatial edges. Before the pulse the 
most convenient free-particle orbitals in (5) are plane 
waves (f>n{x) = ^ (l/_L)e**''"^ where fc„ = 2n-K/L and 
n — —np, ~np - ■ ■ ^np — l,nF with np = {N — l)/2. 
Let (f>n{x, t) be the solution of the TDSE with the Hamil- 
tonian (12) reducing to the above (j)n{x) just before 
the pulse. Then the solutions just after the pulse are 
(/)n{x,0+) = 4'n{x)e'-^^^^^^ . The potential gradients at 
the pulse edges impart momentum kicks to the particles 
there which induce both compressional waves propagat- 
ing at the speed, c, of sound and density dips (gray soli- 
tons) moving at speeds \ v\ < c. The expansion of 0„(a;, t) 
in terms of the unperturbed plane waves is evaluated as 



.ix,t)^l{l + e^^<^)-L^ 



E 



(-l)Vn-2£-i(a;)e-'""-^'-i 
2^+1 



— OO 
t 



(13) 



and the time-dependent density is the sum of the abso- 
lute squares of the lowest N of these. Figure |l| shows nu- 
merical simulations obtained using Eq. (13) for = 51, 
t/re = 0.051, and A0 = tt (solid line), and A9 = 0.57r 
(dashed line) : due to symmetry we show only half of the 
ring —L/2 < x < 0, the phase-shift being imposed at 
X = — i/4. Considering times short compared to the 
echo time means that the corresponding results are not 
very sensitive to the periodic boundary conditions, and 
also therefore apply to a linear geometry. The initial 
density profile is flat with a value poL = 51. For both 
phase-shifts two distinct maxima are seen, which travel 
at close to the speed of sound c, and two distinct min- 
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ima, which are analogous to gray sohtons and travel at 
velocities |u|/c < 1. 




FIG. 1. Scaled density p{x,t)L versus scaled position 
around the ring x/L for = 51, t/r^ — 0.051, and A9 = vr 
(solid line), and A9 — 0.5tv (dashed line). Due to symmetry 
we show only half of the ring — L/2 < a; < 0, the phase-jump 
being imposed ai x = — L / 4. 

In addition, there are also high wavevector oscillations 
which radiate at velocities greater than c. In the case of a 
phase-shift A9 — it, the density is symmetric about x = 
—L/A, whereas for a phase-shift other than a multiple 
of TT the evolution is not symmetric, see the dashed line 
where the global minimum moves to the right in reponse 
to the phase-shift. 



eral features, the generation of gray solitons and density 
waves, agree with those of the GP theory, but here arise 
out of the exact many-body calculation. 

In conclusion, using our exactly-soluble ID model we 
hope to have shown that the dark solitonic features of 
atomic BECs normally described within the mean-field 
GP theory arise naturally from consideration of the ex- 
act linear many-body theory for times less than the echo 
time. An advantage of this approach is that it is number- 
conserving and does not rely on any symmetry-breaking 
approximation. In addition, long time dynamics such as 
collapses and revivals are accounted for Q. A detailed 
comparison between our results and current experiments 
is not possible as they do not conform to the conditions 
for a ID system. However, some estimates are in or- 
der to set the appropriate time scales: If we consider 
^^Rb with a ring of circumference L = 100 /im, and a 
high transverse trapping frequency uj± — 27t x 10^ Hz, 
then we are limited to atom numbers N < 300 so 
these are small condensates. We then obtain — 4.6 
s, and Te ~ 90 ms for N = 51. Finally, we remark that 
since our approach relied on the mapping between the 
strongly-interacting Bose system and a non-interacting 
"spinless Fermi gas" model, this suggests that dark and 
gray solitons should also manifest themselves in the den- 
sity for the ID Fermi system. Although real fermions 
have spin, the interactions used here to generate solitons 
were spin-independent. 

This work was supported by the Office of Naval Research 
Contract No. N00014-99-1-0806. 




FIG. 2. Dark soliton velocity \v\/c scaled to the speed of 
sound c as a function of phase-shift A6/n for A'' = 51. 

In Fig. H we plot the calculated velocity of the global 
density minimum relative to the speed of sound for a va- 
riety of phase-shifts A9. The basic trend is that larger 
phase-shift means lower velocity, in qualitative agreement 
with recent experiments ||^,^ , but there is a sharp veloc- 
ity peak at A9 ~ O.SStt: This peak results from the cross- 
over between two local minima in the density. These gen- 
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